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From an analysis of the specific heat of one of the cuprate superconductors it is shown, that even 
if a large part of the experimental specific heat associated with the superconducting phase transition 
is due to fluctuations, this part must be counted when one tries to extract the condensation energy, 
Econd, from the data. Previous work where the fluctuation part was subtracted, has resulted in an 
incorrect estimation of Econd- 
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In conventional metals superconductivity has been un- 
derstood to result from an effective attractive interaction 
between electrons. This simultaneously causes a reduc- 
tion of the interaction energy and an increase of the ki- 
netic energy when the material becomes superconduct- 
ing. For high Tc superconductors, and in particular for 
the cuprate family, it has been proposed that the op- 
posite situation could exist, where the superconducting 
state is accompagnied by a reduction of the charge car- 
rier kinetic energy ^, |[ ^. Experimental investiga- 
tions of the kinetic energy component perpendicular to 
the superconducting planes of the cuprate high Tc su- 
perconductors have previously demonstrated that the ki- 
netic energy reduction perpendicular to the layers is far 
too small to account for the condensation energy ||, ^, 
which ruled out interlayer tunneling (ILT) as a mech- 
anism of superconductivity|, |, |, fll O . Al- 
though evidence was subsequently reported||13|, |14| that 
the c-axis kinetic energy is reduced in the superconduct- 
ing state in a number of cases, the amount of energy is 
orders of magnitude smaller than earlier estimates of the 
condensation energy]^, ||, |^. 

Later experiments have identified two contributions to 
the internal energy of cuprate superconductors: (i) From 
a re-analyses of inelastic neutron scattering data it was 
concluded that the a6-plane spin-correlation energy was 
lowered by an amount which may be sufficient to account 
for the condensation energy [p5|, and (ii) an even larger 
lowering of a6-plane kinetic energy was measured with 
optical spectroscopy . Since spin-correlations result 
from exchange interactions, which in turn reflect the spin- 
dependent virtual motion of electrons, these two channels 
of condensation energy may have at least in part the same 
microscopic origin. 

However, the analysis providing the condensation en- 
ergy from specific data has been questioned in 1999 by 
Chakravarty, Kee, and Abrahams[0 (hereafter CKA), 
who stated that "the attempt to extract the condensa- 
tion energy from the specific heat data runs into ambi- 
guity, except within a mean field treatment. In the pres- 
ence of fluctuations, superconducting correlations, which 
can primarily be of in-plane origin, contribute to the en- 



ergy and significantly to the specific heat of the normal 
state." In order to resolve this ambiguity, CKA pro- 
posed to "subtract the fluctuation effects and to use the 
remainder as an effective specific heat from which to ex- 
tract the c-axis contribution to the condensation energy. 
The rationale is that free energy can be decomposed into 
a singular and a non-singular part. The universal sin- 
gular part is more sensitive to collective long-wavelength 
fluctuations, while the non-singular part is dominated by 
short distance microscopic pairing correlations." Sub- 
tracting the fluctuation led CKA to a value of the con- 
densation energy which was 40 times smaller than the 
value obtained in Refs. ^, ||, ||. 

Here we demonstrate that the analysis of CKA was 
internally inconsistent. If carried out correctly, the sub- 
traction of fluctuation energy makes only a factor of 2 
to 3 difference compared to Refs. ^, ^, ^. Moreover, we 
show that it is overwhelmingly natural to count also the 
fluctuation contribution in the condensation energy. 

The condensation energy, Econd is the internal energy 
of the equilibrium phase relative to the internal energy 
of the normal state. The former is the experimentally 
observed phase, which is superconducting for T < Tc, 
whereas the latter corresponds to the state where all su- 
perconducting correlations have been suppressed in the 
sense that the two-particle correlation function tends to 
zero as a function of the 'center-of-mass' variable over 
a range no greater than a few times the inter-particle 
spacing. In the remainder of this paper we will use the 
sub-index n to indicate the thermodynamic quantities 
corresponding to this normal state. In any supercon- 
ductor long range phase-coherence is only present for 
T < Tc- In BCS theory long range phase coherence 
and pair-correlations become non-zero simultaneously for 
T < Tc- Knowledge of equilibrium and normal speciflc 
heat for T < Tc then suffices to determine Econd- On 
the other hand, pair-correlations can in principle still ex- 
ist for temperatures above the transition temperature, 
and indeed such correlations are often associated with 
the pseudo-gap phenomenon in underdoped cuprates. A 
measurement of the internal energy, released when the 
superconducting state is formed out of the normal state. 
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FIG. 1: Open circles: Electronic specific heat data of 
T12201. Dashed curve: Singular (fluctuation) contribution, 
parametrized as Csing = 5±A, where t = \1 ~ T /Tc\ {Tc — 
79 K,g+ = 2.38, <?_ — 0.74). Solid curve: Electronic specific 
heat with the singular contribution subtracted. 

should now also include the pair-correlations which al- 
ready exist above the superconducting phase transition. 
Since our discussion is most easily formulated in terms 
of the entropy, let us remind the reader that the entropy 
follows uniquely from the specific heat according to the 
relation 

siT)^£^dr (1) 

// the temperature dependence of the specific heat is 
known in equilibrium and in the normal state, the free 
and internal energy differences can be calculated directly 
using the relations 

/>oo 

F„(r)-F(T)= / {Sn{T')-S{T'))dT' (2) 

JT 
/•oc 

Er,{T)-E{T)^ / {C{T')^C,,{T'))dT' (3) 
Jt 

The integration limits ensure that Iuht^oo FiT) ~ 
FniT) and Ivcht^oo E{T) ~ £'„(T), in accordance with 
the assumption that for T ^ c» all superconducting cor- 
relations vanish. The condensation energy corresponds to 
the zero-temperature energy difference Econd = En{0) — 
E{0) — F„(0) — F{0). This positive energy can be ob- 
tained either by integrating the specific heat difference, 

/•OO 

Econd ^ / {C{T)-Cn{T))dT (4) 

Jo 

or by integrating the entropy difference 

/>OC 

Econd = / {Sn{T)-S{T))dT (5) 

Jo 

The fact that Eq. | and Eq. | should give the same 
value provides, as we will see, an important consistency 
check. CKA[Q questioned the analysis providing the con- 
densation energy, and provided a different analysis where 
a Gaussion fluctuation contribution was subtracted from 
the experimental data. In Fig. |l| the fit obtained by CKA 



to the specific heat of T12201 to 2D Gaussian fluctuation 
plus non-singular parts is reproduced. We will indicate 
with an asterisk thermodynamic quantities from which 
the fluctuations have been subtracted. CKA calculated 
both results with (£;„(0) - E*{0) = 25 mJ/g-at.) and 
without (£'„(0) - E{0) = 825 mJ/g-at.) subtracting the 
singular part. The latter corresponds to the earlier esti- 
mates in Refs. ^, ^, ^. Following the basic idea of the 
ILT theory, CKA now equate the 'subtracted' value of the 
condenstaion energy to the decrease of the c-axis kinetic 
energy in the superconducting state, SK. Finally they 
use the standard relation |, |, |, |, |l^, |ll|, H H |l|] 
between SK and the c-axis penetration depth Ac, 

8A! " 2ld??'^^' 

Their resulting estimate of Ac is much larger than the 
value Ac — 1/im derived previously in the context of ILT. 
Vice versa the implication was, that the condensation 
energy is approximately 40 times smaller than earlier es- 
timates, which had not corrected the speciflc heat for the 
singular part. 

The main source of this huge difference is the diffi- 
culty in determining the 'normal' thermodynamic quan- 
tities. The experimental questions include not only how 
but also whether the 'normal' state can be reached as 
a function of temperature, or indeed anything else. In 
other words, whether or not a parameter exists that can 
be tuned to lower the energy of the normal state below 
that of the equilibrium superconducting state so that it 
can be accessed without losing its fundamental character. 
Quite generally such a parameter does not always have 
to exist. For a weak coupling superconductor a mag- 
netic field would suffice. In most cuprate superconduc- 
tors the required magnetic fields are beyond experimental 
reach, but Zn ions substituted for planar Cu may serve as 
an alternative for suppressing superconductivity jl|, p^ . 
However, for the cuprates there is reason to believe that 
several 'normal' states are competing with the supercon- 
ducting one (eg stripe, DDW, normal). In this case the 
field (or Zn-doping) required to mute superconductivity 
could be enough to rearrange the order between these 
'normal' states, thus revealing the 'wrong' one when su- 
perconductivity gets suppressed. Thus we are confronted 
with the difficult situation, that the normal state entropy 
is not an experimental quantity, and can only be deter- 
mined based on theoretical considerations and/or by ex- 
trapolating the normal state dependence as was indeed 
done in Ref. 0, providing us, as we have seen, with es- 
timates of condensation energy differing by a factor of 
40. 

However, the situation isn't as bad as it looks. We 
can use our knowledge of thermodynamics to constrain 
the possible behaviour of S'„(T): Both S{T) and 5'„(T) 
are subject to the 2d and 3d law of thermodynamics. We 
can also use the reasonable assumption, that for tempera- 
tures high enough, all superconducting correlations cease 
to exist, causing S'„ and S to become equal in that limit. 
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FIG. 2: Experimental entropy versus temperature. Dashed 
curve: Normal entropy with 7 = 0.576 mJ/G-at.K^ 

The circumstance that F{T) corresponds to the equilib- 
rium state, impUcs that for any temperature Fn(T) has 
to be larger than F{T). The corresponding constraints 
on the entropy are, in the same order, 

dSn/dT > 
Sn{0) = S{0) = 
Snioo) = S{oo) 

Sn{T')dT' < S{T')dT' + E,ond 

In Fig. ^ the entropy is plotted as a function of temper- 
ature. The condensation energy in this plot corresponds 
to the area between the equilibrium entropy [S{T): open 
circles) and the normal entropy (S'„(T): dashed curve). 
A condensation energy which is 40 times smaller than 
the earlier estimates of the condensation energy^, ||, |j, 
would require that the area enclosed by 5„(r) and S{T) 
is also 40 times smaller than in Fig. ||. Due the con- 
straints of Eq. ^ Sn{T) would almost coincide with 
S{T). Sn(T) then has a sharp bend at Tc, which would 
correspond to a pathological phase transition of the nor- 
mal state in close proximity to the real superconducting 
phase transition. 

Let us now return to the analysis by CKA. If we 
calculate the entropy by integrating C*{T)/T with the 
singular part subtracted, we obtain the curve indicated 
with open circles in Fig. ^. Above the phase transi- 
tion S*{T) approaches a linear temperature dependence 
with a negative offset, given by S*{T) ^ - So. We 
notice, that the corrected entropy S*{T) and the nor- 
mal entropy do not merge above Tc, causing the integral 
E*cond = ir i^n{T) - S*{T)) dT to diverge. We see, that 
subtracting the singular part has the effect that two pro- 
cedures by which the condensation energy can be calcu- 
lated provide opposite results: 

Econd = Jo" {C{T) - Cn{T)) dT 25mJ/g-at. 

Econd = iT iSn{T) - S{T)) dT^^ 

This inconsistency does not arise in the earlier esti- 
mates of the condensation energy]^, ^, ||]: The fact that 
Sn{T) and S{T) merge above T^ in Fig. || guarantees 
that the same value of Econd is obtained with the two 
different formulas. 

Following a recent suggestion by Abrahams the 
CKA-analysis can be improved by adopting a specific 



FIG. 3: Open circles: Corrected entropy versus temper- 
ature S*(T) = S{T) - SsingiT), following the procedure 
of Ref. [ij where the singular contribution was subtracted. 
Dashed curve: Normal entropy, S„{T), same as in Fig.^. Dot- 
ted curve: 'Normal' entropy, S*{T) fitted to the T > Tc- 
region of the corrected entropy curve. Solid curve: Singular 
contribution to the entropy SsingiT). 

heat exponent rj different from 1 below Tc- We must 
also take into account that the entropy above for T > Tc 
approaches S*{T) = 7T — Sq. A simple analytical ex- 
pression with these two properties is 

^:(r) = 5o|(^l+(^^y) ^'-ij (8) 

We have tried to fit different values for rj. If we choose 
rj > 2.5 the 'normal' entropy S*{T) becomes smaller than 
the experimentally measured S(T) for temperatures be- 
low 40 K. Although such a possibility can not be excluded 
by the thermodynamical constraints of Eq. |7[ from a 
microscopic perspective it looks suspicious that the en- 
tropy of the (gapped) superconducting state could exceed 
that of the normal state. To avoid this, we have adopted 
the value 77 = 2.0. The best fit in the region between 
Tc and IIQK was obtained with Sq = 34.85 mJ/g-at.K 
and To = 50K. In Fig.^ we display the singular con- 
tribution to the entropy SsingiT), the corrected entropy 
S*iT) = SiT) — SsingiT), and the normal state entropy 
S'*(T). Note that if indeed it would be justified to sub- 
tract the singular contribution, the presence of a nega- 
tive offset in the entropy implies that this 'normal' state 
would have a pseudo-gap fs], pj\ . 

In Fig. I we display both S'„(T) - SiT), and the quan- 
tity which corresponds to the improved version of the 
CKA-analysis, S*iT) - S*iT). We see, that with the im- 
proved CKA-analysis the difference entropy between the 
normal and superconducting state indeed becomes zero 
above Tc- We carried the analysis of the specific heat 
one step further, and calculated the free energy differ- 
ence curves. These are displayed in Fig. |^. The cor- 
rected CKA-analysis thus provides a condensation en- 
ergy of 0.35 J/g-at. This value is about 30 percent of 
the condensation energy from direct integration of the 
total entropy of Fig.^ which gives Econd = 1-3 J/g-at.. 
We see, that even if it would be justified to subtract the 
fluctuations, the correct analysis would still reduce the 
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FIG. 4: Open circles: Normal state entropy minus experimen- 
tal entropy. Solid curve: Singular contribution to the entropy. 
Dashed curve: 'Normal' entropy Sn{T), minus the corrected 
entropy S*(T). 
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FIG. 5: Integrated entropy differences of Fig. ^ Open cir- 
cles: JgiSniT') - S{T'))dT'. Dashed curve: j^f (S;(r') - 
S* {T'))dT' . Solid curve: Singular contribution to the free 
energy. 

condensation by a factor of 2 or 3, not by a factor of 40 
as stated in Ref. |]. 

However, removing the singular contribution from the 
experimental data comes with a penalty: As a reminder 



that a fluctuation contribution has been subtracted in 
the CKA-analysis, we have displayed in Figs. ||, || and ^ 
its contribution to the entropy and the free energy. We 
see, that the fluctuation entropy has a conspicious step at 
Tc, which indicates that these fluctuations are intimately 
connected to the superconducting phase transition. The 
fact that the fluctuation speciflc heat is strongly peaked 
at the phase transition, implies that there is an additional 
reduction of the internal energy due to the fluctuation 
contribution. Hence it would seem to be overwhelmingly 
natural to count this part when one estimates the energy 
by which the superconducting state is stabilized relative 
to competing (non-superconducting) phases. 

In conclusion, we have shown, that in the original anal- 
ysis of CKA, i.e. subtraction of a fluctuation contribu- 
tion to the specific heat, a corrected entropy below Tc 
was used which does not match the value used above 
Tc. This internal inconsistency was the main reason why 
the condensation energy was estimated a factor of 40 
smaller than the value obtained in earlier publications. 
This problem could have been fixed by letting the spe- 
cific heat exponent be less than one at low temperatures. 
We have repeated the CKA-analysis with this fix, result- 
ing in a condensation energy which is ten times as large. 
However, it seems overwhelmingly natural to include the 
contribution of the fluctuations in the analysis, because 
experimentally the fluctuation contribution to the inter- 
nal energy appears to be intimately linked to the super- 
conducting phase transition. This results in a conden- 
sation energy of approximately 1.3 J/g-at. for optimally 
doped T12201. 
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